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O ■ ABSTRACT 

NGC1313 X-2 was among the first ultraluminous X-ray sources discovered, 
and has been a frequent target of X-ray and optical observations. Using the 
HST/ACS multi-band observations, this source is identified with a unique coun- 
terpart within an error circle of &!2. The counterpart is a blue star on the edge of 
a young cluster of < 10 7 years amid a dominant old stellar population. Its spec- 
tral energy distribution is consistent with that for a Z=0.004 star with 8.5 M & 
about 5 x 10 6 years old, or for an 07 V star at solar metallicity. The counterpart 
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exhibited significant variability of Am = 0.153 ± 0.033 mag between two F555W 
observations separated by three months, reminiscent of the ellipsoidal variability 
due to the orbital motion of this ULX binary. 

Subject headings: Galaxy: individual(NGC 1313) — X-rays: binaries 

1. INTRODUCTION 
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A class of luminous X-ray point sources, usually called ultraluminous X-ray sources 
(ULXs), has been established as non-nuclear sources with X-ray luminosities in the range 
10 39 -10 41 erg/sec. ULXs were discovered in nearby galaxies first by EINSTEIN (e.g., Fab- 
biano et al. 1989), ROSAT (e.g., Colbert & Mushotzky 1999), ASCA (e.g., Makishima et al. 
2000) and now the Chandra and XMM X-ray observatories (e.g., Kilgard et al. 2002; Swartz 
et al. 2004). ULXs have stirred great excitement in the high energy community because 
some may be intermediate mass black holes (IMBHs) of 10 2 -10 5 M & (Miller & Colbert 2004). 
This is because their luminosities exceed the Eddington limits for 10 M Q black holes, and 
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IMBHs are required if the sources are emitting isotropically at their sub-Eddington levels. 
Alternatively, ULXs could be stellar mass black holes with beaming effects, either mild geo- 
metric beaming (King et al. 200f) or relativistic beaming (Koerding et al. 200f), or with a 
radiation pressure dominated magnetized accretion disk where the photon-bubble instability 
operates and the radiation can exceed the Eddington limit by a factor of 10-100 (Begelman 
2002). 

NGC1313 X-2, located about 6' south of the nucleus of NGC1313 at a distance of 
3.7 Mpc (Tully 1988), was among the first ULXs discovered by EINSTEIN (Fabbiano & 
Trinchieri 1987). It has been a frequent target of ROSAT, ASCA, Chandra and XMM- 
Newton observations, and has shown clear variability of up to a factor of 2 on a timescale of 
months (c.f. Zampieri et al. 2004 and reference therein; Feng & Kaaret 2006). The ULX has 
a maximum X-ray luminosity of 3 x 10 40 erg/sec in 0.3-10 keV, and has switched between 
high flux state with hard spectrum and low flux state with soft spectrum, unlike normal 
Galactic black hole X-ray binaries (Feng & Kaaret 2006). Miller et al. (2003) discovered, 
along with a power-law component, a cool accretion disk component of kTi n R3 160eV in its 
X-ray spectrum, suggesting an IMBH of ~ 1O 3 M . This IMBH interpretation, however, is 
not unquestionable and the spectrum can also be fitted by a stellar mass black hole with, e.g., 
a hot accretion disk (kT in ~ 2.6keV) plus a power-law component (Stobbart et al. 2006). 

NGC1313 X-2 has been studied extensively in the optical to understand its companion 
and environments. A deep H a image reveals the ULX as located in an elongated 25" x 17" 
supershell, which showed strong [S II] and [O I] lines with a FWHM of 80 km s _1 (Pakull 
& Mirioni 2002). Zampieri et al. (2004) identified the ULX with a point-like counterpart 
on an R image taken with the ESO 3.6m telescope. Their subsequent efforts (Mucciarelli et 
al. 2005) resolved the counterpart into two sources CI and C2, with CI claimed as a BO-09 
main sequence star of 20M Q and C2 as a G supergiant of ~ 10M o . Ramsey et al. (2006) 
studied HST ACS/WFC images and identified the ULX with CI, which they claimed an 
early B giant. Pakull et al. (2006) studied ESO/VLT and SUBARU observations, and found 
the counterpart less than 8 M in an environment of 40-70 Myrs old. They detected the He 
II A4686 emission line in two blue spectra of the counterpart taken 3 weeks apart, with a line 
shift of 300 km/s. The He II A4686 line indicates the high ionization state for He, probably 
due to photoionization by the strong X-ray emission, or photoioniztion by Wolf-Rayet stars 
or shocks (e.g., Pakull & Mirioni 2002). 

In this paper we report our studies of NGC1313 X-2 with the HST ACS observations. 
In section 2, we describe the observations utilized and the analysis procedures to identify 
the counterpart and to compute the photometry. The ULX is identified with a unique 
counterpart within an error circle of 0"2. Its spectral energy distribution is consistent with 
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that for an 07 V star, and it showed variability above 4<r between two F555W observations. 
The counterpart has an age of 5 x 10 6 years, and is located next to a loose cluster of blue 
stars. We discuss the implications of these results in Section 3. 

2. DATA ANALYSIS AND RESULTS 

NGC1313 X-2 was observed with HST ACS in four filters F330W, F435W, F555W 
and F814W on 2003-11-22, and in filter F555W on 2004-02-22 for program GO-9796. All 
observations (Table 1) were downloaded from the HST archive, and calibrated on the fly 
with the best calibration files. 



2.1. Astrometry 

The drizzled ACS/WFC images were used to identify the ULX. These images were gen- 
erated by co-adding multiple exposures with the same filter and pointing with the DRIZZLE 
package (Fruchter & Hook 2002), which was able to reject most cosmic rays plaguing single 
exposures. The geometric distortion of the camera was corrected to better than 0.01 pixels 
(ACS handbook). The X-ray positions were taken from a 15 ksec Chandra ACIS-I observa- 
tion (ID 3550), which was downloaded from the Chandra Data Archive and processed with 
CIAO 3.3.0. WAVDETECT (Freeman et al. 2002) was run to detect point sources, and two 
other X-ray sources (X3 and X4) were detected within 1' of the ULX (Table 2). 

The direct registration of Chandra images on HST images yields a nominal error circle 
of ~ 0'/8, as a quadrature sum of the absolute position uncertainties from Chandra images 
(0"6, Aldcroft et al. 2000) and HST images (0"5, HST data handbook). Luckily, the ULX 
position on WFC images can be greatly improved with the help of X4, which was identified 
with the nucleus of a background galaxy with a size of 0'/5 as in Figure la. The correctional 
shift between the Chandra and HST positions of X4 was applied to the Chandra position 
of the ULX to obtain its HST position (Table 2). The positional error from this procedure 
included the ACIS plate scale variation (~ 0'/03 for the 36" separation between X4 and the 
ULX), the centroiding errors for X-ray sources (0'/17 for X4 and O'/Ol for the ULX; Table 2). 
These amounted to an error circle of < 0'.'2 around the ULX as shown in Figure lb. Within 
this error circle was only one point-like object, which was the counterpart candidate CI in 
Mucciarelli et al. (2005), and another candidate C2 was unambiguously excluded. While 
the highly absorbed X-ray spectrum and highly variable light curve of X4 were consistent 
with its identification as the nucleus of the background galaxy, it nevertheless could be any 
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source in the dusty spiral arms. In this situation, the positional error had contribution from 
the size of the background galaxy, and could be as large as Of! 5. 

X3 had a very soft spectrum and was identified with a foreground star that was highly 
saturated on the WFC images. The optical position of X3 could be estimated roughly as the 
crossing of the two spikes as listed in Table 2. Despite the uncertainty (~ Of! 5) of this optical 
position, the ULX position obtained through X3 was consistent with that obtained through 
X4, as shown on Figure lb. Note that the HST coordinates were offset ~ 2'.'3 relative to the 
Chandra coordinates. Such an offset was larger than the nominal error of 0f!8, but was not 
uncommon, and was real given the successful identifications of both X4 and X3. In fact, in 
the study of another ULX in M81, the HST and Chandra positions for a conspicuous source, 
SN1993J, were offset by ~ 2" (Liu et al. 2002). 



2.2. Photometry 

Photometry for the ULX counterpart and nearby objects was computed on the drizzled 
images with the IRAF/DAOPHOT package. The drizzled images, originally given in unit 
of e~/s, were first converted to e~ per pixel by multiplying the total exposure times (Table 
1). An aperture of 3 pixels was used to extract the flux in the task PHOT, and the aperture 
corrections as taken from Sirianni et al. (2005) were applied to obtain the magnitudes for 
an "infinite" aperture. A fitting radius of 3 pixels was adopted to fit PSFs to objects. 
We computed the magnitudes both in the VEGAmag and the STmag systems, with the 
zeropoints taken from Sirianni et al. (2005). The photometry results from DAOPHOT 
were compared to those obtained with the HSTPHOT package (Dolphin 2000), which were 
consistent with each other to better than ±0.1 mag. 

Variability of the counterpart was studied with the two F555W observations in 2003 and 
in 2004, which were separated by three months. The sources from the two observations were 
cross identified, which led to 399 bright sources within 10" of the counterpart. As shown in 
Figure 2, most sources did not exhibit significant variability between the two observations, 
confirming the stability of our photometry procedures. For 318 (i.e., 80%) of these objects, 
the variations were less than 2o"A m , with o^ m = ^5m\ + 5m\. Twelve sources exhibited 
variability above the 4o"A m significance level, and the counterpart was one of them with 
ArriF555w = 0.153 ± 0.033, an amplitude of 15% at a significance of 4.6o"Am- 

The multi-band observations in 2003 were used to study the spectral energy distribution 
of the ULX counterpart. The flux densities at the central wavelengths for F330W, F435W, 
F555W and F814W were calculated from the STmag in these bands, and plotted in Figure 
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3 with corrections for different extinction values based on the extinction law prescribed 
by Cardelli, Clayton & Mathis (1989). The observations were compared to the absolute 
magnitudes and colors for blue stars of 05V, 07V, 09V and B1III spectral types as taken 
from Schmidt-Kaler (1982). If the counterpart is only corrected for the Galactic extinction 
of 0.11 mag toward NGC1313 (Schelegel et al. 1998), its F435W and F555W magnitudes 
are roughly consistent with those for an 09V star, but its spectral slope is shallower and 
its F814W magnitude is much brighter. If the counterpart is corrected for an extinction 
of E(B-V) = 0.33 mag, its overall spectral shape matches perfectly that for an 07V star 
[My = —5.2 mag, ~ 3OM , ~ 9R Q ). The spectral shape will be over corrected if we apply 
E(B-V) = 0.44 mag, i.e., = 2.7 x 10 21 cm -2 as inferred from the power-law fit to the 
X-ray spectrum (Miller et al. 2003). The Bl III star has B and V magnitudes consistent 
with the counterpart, but its U magnitude is lower than the observation as shown in Figure 
3. Note the properties for the MK spectral types are typical for stars with solar metallicity 
in our Solar neighborhood, therefore different from low metallicity stellar models. 



2.3. Color-Magnitude Diagrams 

The color-magnitude diagrams for the ULX and nearby bright stars were constructed 
to estimate their ages. Following the advice of Sirianni et al. (2005) against transforming 
ACS VEGAmag to standard magnitudes, we constructed the diagrams with VEGAmag and 
compared to isochrones computed in the HST ACS/WFC VEGAmag syste for the Padova 
stellar models (Girardi et al. 2002). Studies of the interstellar medium of NGC1313 indicate 
a rather low metallicity of about 0.1-0.2 Z Q (Ryder 1993), thus isochrones for Z=0.004 were 
used to overlay on the diagrams in Figure 4. The magnitudes and colors in these diagrams 
were corrected for the Galactic extinction of E(B-V)=0.11 mag toward NGC1313. 

Two stellar populations were clearly seen from the color-magnitude diagrams, a young 
population with ages less than a few 10 7 years and an old population with ages 3 — 30 x 10 s 
years. The ULX counterpart belongs to the young stellar population. By comparing its 
position on the diagrams with those of the stellar models, we found that the counterpart 
has an age of 10 7 years from its F435W and F555W magnitudes after corrections for E(B- 
V)=0.11 mag, but has an age of 3 x 10 7 years from its F555W and F814W magnitudes. The 
two ages have smaller discrepancies for larger extinction corrections, and converge for E(B- 
V) = 0.33 mag, for which the counterpart has an age of 5 x 10 6 years, an initial/current mass 
of ~ 52/8. 5M , and a radius of ~ 7R & . If this extinction is applied to all stars, the young 
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population has ages of 5 — 10 x 10 6 years, and the old population has ages of 2 — 20 x 10 8 
years. 

The two stellar populations are obvious on the accompanying 50" x 50" true-color image 
centered at the ULX counterpart (Figure 7). The image was constructed from the F435W, 
F555W, and F814W observations in 2003, and the bright cosmic rays not removed by the 
DRIZZLE package were manually removed. The field was dominated by old red faint stars, 
with scattered young blue stars. A few tens of bright blue stars were loosely clustered within 
an irregular region of 10" x 5" (i.e., 180pc x 90pc) as outlined on the image. The counterpart 
is the marked bright blue star located to the east-south edge of this loose cluster. The cluster 
region is much larger than for typical Galactic open clusters (< 10pc), and the member stars 
may be the products of one star formation episode that are not bounded gravitationally. 
We note that this loose cluster consists of only blue stars, and no red supergiants. For the 
adopted stellar models with Z = 0.2Z Q , only stars less massive than 2OM can evolve to 
red supergiants in ~ 10 7 years or longer. The absence of red supergiants indicates the open 
cluster is younger than ~ 10 7 years, regardless of the extinction correction. This is consistent 
with the results from comparison with isochrones. 



3. DISCUSSION 

Different efforts to identify the counterpart for NGC1313 X-2 have now converged upon 
the same blue star. By registering X-2 onto an ESO 3.6m R image with respect to SN1978K 
at ~ 5' away, Zampieri et al. (2004) identified a point-like object within the 1"4 error circle, 
which was later resolved into two stars, CI and C2, on ESO VLT images (Mucciarelli et al. 
2005). Ramsey et al. (2006) used stars in the USNO B1.0 catalog to determine astrometric 
solutions for the HST ACS/WFC images to an rms accuracy of / /5, and identified CI as 
the counterpart. In this paper, we registered X-2 onto the HST ACS/WFC images with 
respect to two nearby X-ray sources that were both identified optically, and confirmed CI 
as the counterpart within an error circle of 0"2. The chance for this blue star to be a bright 
interloper is negligible (< 10~ 2 ) given the small size of the error circle and the local stellar 
surface density. Furthermore, the blue star exhibited significant variability between two 
F555W observations, as expected if this blue star is the secondary of the X-ray binary. We 
conclude that the blue star CI is the true counterpart for NGC1313 X-2. 

Ramsey et al. (2006) studied the counterpart with the transformed B and V magnitudes 
and claimed it a 7-9 M Q early B giant. They studied the nearby stars without any extinction 
correction and concluded that the stellar population was at least 10 7 years old without any 
O stars. In this work, two more bands, F330W and F814W, were included in the spectral 
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energy distribution of the counterpart, and it was found that an early B giant such as B1III 
does not match the overall spectral shape of the counterpart. The counterpart is brighter 
and hotter than a B1III giant, and should have a spectral form of v~ 2 in the observed optical 
bands, i.e., the Jeans tail for a blackbody spectrum. To correct the observed spectral shape to 
z/~ 2 , an extinction value of E(B-V) = 0.33 mag is required. This value is, reasonably, lower 
than the value derived from X-ray absorption (0.44 mag); it is higher than the Galactic 
value (0.11 mag), suggestive of additional extinction in NGC1313 and the very vicinity of 
NGC1313 X-2. With this extinction correction, the counterpart has magnitudes and colors 
for a, Z = 0.2Z & star with an age of 5 x 10 6 years, a current mass of ~ 8.5M , and a radius 
of ~ 7Rq. Incidentally, the corrected magnitudes and colors are also consistent with those 
for an 07V star with Z = Z Q , a mass of 30 M Q , and a radius of 9 R & . This identification of 
the ULX as an 07V star joins previous ULX identifications as massive OB stars (e.g., Liu et 
al. 2002; Liu et al. 2004; Kaaret et al. 2004) to corroborate the ULX scenario in which the 
massive secondary is overflowing its Roche Lobe to supply the high accretion rate required 
for the high X-ray luminosity (king et al. 2001). 

Pakull et al. (2006) speculated that the optical light from the X-ray illuminated the 
accretion disk dominated that from the secondary which they claimed a B giant. They 
argued that the X-ray illuminated accretion disk can reach My ~ —5 mag for a stellar mass 
black hole with Lx/Le ~ 100 and P ~ 200 hours based on the van Paradijs & McClintock 
S — My diagram (1994) derived for Galactic low mass X-ray binaries, if the relation still 
holds for Lx > 10 40 erg/s. For an X-ray illuminated accretion disk, the temperature depends 
on the radius R roughly as i? -a5 (e.g., Frank, King & Raine 2002; Raymond 1993), and 
decreases to about 4000 K at R < 9R@. The resulted spectral energy distribution takes 
a form of v~ x in the observed bands, redder than the O star spectrum (z/ -2 ) by a color 
difference A(B — V) ~ 0.25 mag^l The presence of significant accretion disk light would 
modify the classification of the secondary and its estimated properties. In Section 2, we 
corrected the observed spectral energy distribution to that for an 07V star (My = —5.2 
mag) by an extinction of E(B-V) = 0.33 mag, an extra E(B-V)=0.22 mag in addition to the 
Galactic E(B-V)=0.11 mag. Alternatively, the observed spectral energy may be interpreted 
as 12% from the secondary {v~ 2 ) plus 88% from the accretion disk (z/" 1 ) corrected only by 
the Galactic extinction E(B-V)=0.11 mag, or 56% from the secondary plus 44% from the 
accretion disk corrected by an extinction of E(B-V)=0.22 mag. In the first alternative, the 
secondary would have My = —2.9 mag, similar to a B1V star (My = —3.2 mag) or slightly 



2 For the O star spectrum, B — V = 2.5lg(f v /f B ) - 2.5lg(fvo/ fm) = 2.5lg(X v /X%) - 2.5lg(f VO / '/bo); 
for the accretion disk spectrum, B — V = 2.5lg(fv / f B ) — 2.5lg(fvo/ /bo) = 2.blg(Xv / X B ) — 2.5lg(fvo/ /bo)- 
Thus A(B - V) = 2.Ug{X v /X B ) ~ 0.25 mag for Ay = 0.55 nm and X B = 0.435 nm. 
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later. In the latter, the secondary would have My = —4.6 mag, similar to a B1III giant 
(My = -4.4 mag). 

We suggest the accretion disk light does not dominate the secondary light. The first 
of the above alternatives is not reasonable because there must be some extinction intrinsic 
to NGC1313. The second alternative is not true either, or else the composite spectrum 
would show a drop in F330W due to the cool B1III giant, which we did not observe. In 
fact, the same S — My relation predicts a My ~ —2.2 mag for the accretion disk around an 
IMBH of rsj 10 2 M Q and a period of ~ 50 hours as shown later, < 10% of the total observed 
light. Copperwheat et al (2006) carried out detailed simulations for seven ULX systems, 
and found that the light from the disk, depending on factors such as the disk size and the 
inclination angle, is generally much less than that from the secondary. Similarly, we expect 
that the optical light from the accretion disk is insignificant as compared to the massive 
secondary star. Nonetheless, the relative fraction of the optical light from the accretion disk 
can be measured with multiple simultaneous optical and X-ray observations. Assume the 
total optical light L opt = L secon dary + Ldisk, with L secan ^ ary constant in all observations, and 
Ldisk = C x Lx in linear proportion to the X-ray luminosity Lx- In principle we can solve or 
best fit for L seconi i a r y and the linear factor C if we have multiple (2+) observations in which 
the X-ray luminosities vary significantly. 



3.1. The association with a young cluster 

We emphasize that there is a young stellar population of < 10 7 years in additional to 
the apparent old population around NGC1313 X-2. Red stars from the old population are 
dominant in this sky region, while blue stars are rarer. Careful inspection of the accom- 
panying true-color image shows that these blue stars are loosely clustered in a ~ 10" x 5" 
region that is not gravitationally bounded. The member stars may have been generated from 
the same star formation episode that is still active. Given the low surface density of this 
young cluster, it is easily missed, leading to claims that the whole stellar population is old. 
This is particularly true for short exposures and low resolution images, and may lead to an 
underestimate of the ULXs that are identified with young massive stars. 

There are implications on the formation of NGC1313 X-2 given its location nearby, but 
not in the center of, a young cluster that is not gravitationally bounded. If NGC1313 X-2 has 
radiated below its Eddington luminosity, its peak luminosity of 3 x 10 40 erg/sec will suggest 
a black hole of > 200M Q . Such an IMBH cannot be formed from massive stars with large 
mass loss due to stellar winds incurred by metal lines, but can be formed through mergers 
in dense clusters. One scenario is the merging of stellar mass black holes in the dense core of 
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globular clusters (Miller et al. 2002). The chance for thus formed IMBHs below 300 M to 
be kicked out of the cluster centers is not negligible, but the time scale is much longer than 
10 8 years to form a 200 M IMBH (Gultekin et al. 2004). Another scenario is the merging 
of massive stars in super star cluster of ~ 1O 6 M (e.g., Portegies Zwart & McMillan 2002). 
Frequent stellar captures and mergers lead to super stars of ~ 10 3 M Q that subsequently 
form IMBHs via direct collapses or supernova explosions in a few million years. The resulted 
IMBHs, however, tend to sit in the centers of super clusters, inconsistent with X-2's location 
outside a young cluster. 

Another formation scenario for NGC1313 X-2 is the merging of proto stars in proto 
clusters as suggested by Soria (2005). Stellar captures and mergers are enhanced by proto 
stellar disks/envelops and low stellar velocities, and super stars of ~ 1O 3 M can form in 
less than a million years in small proto clusters of ~ 10 4 — 1O 5 M . Subsequent explosions 
of these super stars will unbound the proto clusters, and form IMBHs of a few hundred 
M & in weakly bounded or unbounded clusters, usually off-center as found for NGC1313 X- 
2. This formation mechanism works best in low Z environments and in colliding or tidally 
interacting systems as exemplified by the Antennae (Fabbiano et al. 2002), Cartwheel (Gao 
et al. 2003), NGC7714/7715 (Smith et al. 2004) and M99 (Soria et al. 2005). NGC1313 
X-2 is, incidentally, in an environment of low metallicity (~ 0.1 — 0.2Z Q ) in the southern 
satellite regions of NGC1313, which could be a tidally disrupted companion galaxy (Sandage 
& Brucato 1979) or the result of collision of huge HI clouds with the disk of NGC1313 
(Marcelin & Gondoin 1983). 



3.2. Optical variations 

The counterpart exhibited variability of Am = 0.153 ± 0.047 mag (3.3<x) between two 
ACS/WFC F555W observations separated by three months. Most probably, this variability 
reflects the orbital motion of an inclined binary system for NGC1313 X-2. Indeed, the 
variability is about the level of ellipsoidal variability expected for a high-mass secondary 
that approximately fills its Roche lobe. For example, the amplitude is 0.17 mag for an 
05V secondary, a 10 M black hole, L x = 10 40 erg/s, and cosi = 0.5 (Copperwheat et al. 
2005). 

It is unlikely that this variability is caused by the intrinsic variability of a pulsating star 
as the secondary. The observed variability and the blue color may be possible for a (3 Cephei 
star that varies by 0.01 - 0.3 mag with periods of 0.1 - 0.6 days, or a Wolf-Rayet star or a 
Luminous Blue Variable (LBV) which exhibit day-to-day microvariations of 0.1-0.2 mag and 
"normal" variations of 1-2 mag on timescales longer than months. However, the secondary 
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cannot be a Wolf-Rayet star because no prominent emission lines for Wolf-Rayet star were 
shown in VLT observations (Pakull et al. 2006). The secondary cannot be a (3 Cephei star 
or a LBV because the observed luminosity (My ~ —5) is high for a typical j3 Cephei star 
(My - -2), but low for a typical LBV (My 7). 

On the other hand, the optical variability might have sizable contribution from the ac- 
cretion disk in presence of substantial X-ray variability (e.g., flares). Fortunately, the change 
of the optical spectra will enable us to quantify this contribution, because the spectrum of 
an O type secondary is bluer than the accretion disk by a color difference of A(B — V) = 0.25 
mag. For example, if the accretion disk light increases from 5% to 20% of the secondary 
light, the spectrum will get redder by A(B — V) ~ 0.04 mag. Such a color difference can be 
detected by HST observations with 1% photometric accuracy. 

A well designed photometric monitoring program can sample the ellipsoidal light curve 
and obtain the orbital period. The period of NGC1313 X-2 can be calculated as a function 
of the black hole mass (Mi) for the known secondary mass (M 2 ) and size (i^)- The Roche 
lobe size R cr of the secondary is solely determined by the orbital separation a and the mass 
ratio q = M 2 /Mi as R cr = af(q). If the secondary is filling its Roche lobe, i.e., R 2 = R cr , we 
can compute the orbital separation as a = R 2 / f(q), and subsequently calculate the period as 
P 2 = 2n G ^ M ® +M ^ . The predicted periods for different primary masses are shown in Figure 
5, in which f(q) values as tabulated in Kopal (1959) are used. The analytic approximations 
to these tabulated values by Paczynski (1971) and by Eggleton (1983) were also tried, which 
gave periods different from Kopal of less than 5% as shown in Figure 5. 

The predicted periods are largely determined by the secondary properties. Adopting 
Paczynski's approximation for q < 0.8, f(q) = 0.46224(y^) 1 / 3 , the period can be expressed 
as P = a/110/p hr independent of the primary mass, in which p is the mean density of the 
secondary in unit of g cm -3 (Frank, King & Raine 2002). For massive black holes (q < 0.8), 
the period is ~ 44 hours for a secondary of an 07V star (3OM ,9-R ), ~ 115 hours for a 
BOO giant (2OM , 15i? ), and ~ 290 hours for a B0I supergiant (25M , 3Oi? ). Thus, the 
observed period can be used to test possible classifications of the secondary. On the other 
hand, if the secondary mass and size can be determined by other means, the observed period 
can be used to constrain the black hole mass. In the case of NGC1313 X-2, the secondary 
mass and size are determined to be 8.5 M and 7 R & by comparing to Z = O.O2Z stellar 
models. We expect a period of ~ 56 hours for massive black holes, and shorter periods for 
primaries less massive than ~ 15M , i.e., a stellar mass black hole for this ULX (Figure 5). 

The orbital velocities around the mass center for the Roche lobe filling system NGC1313 
X-2 are shown in Figure 6 for different primary masses. For the secondary of 8.5 M and 7 
Rq, the velocity increases from a few tens to a few thousands km/sec for black hole masses 
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from a few to a few thousand M Q . For an IMBH of a few hundred M , the velocity change 
is about 10 3 km/sec, easily detectable with low spectral resolutions of a few A at optical 
wavelengths. The orbital velocity for the black hole (i.e., the accretion disk) decreases from 
~ 300 to a few km/sec for increasing black hole masses. Note the radial velocity is smaller 
than the orbital velocity by sini for an inclination angle i. 

Pakull et al. (2006) detected broad He II A4686 emission lines with FWHMs of ~ 
600 km/s in two ESO VLT spectra of NGC1313 X-2, which they attributed to the X-ray 
illuminated accretion disk. The He II line appeared to have moved by 300 km/sec, and 
may indicate the system has a stellar mass black hole since the orbital velocity change is 
very small for very massive black holes (Figure 6). However, the He II line from an X-ray 
illuminated accretion disk is usually broader, e.g., about ~ 1200 km/s for an inclination 
angle of ~ 30° (Raymond 1993). The emission line will be narrower for smaller inclination 
angles (i.e, face-on), but the full amplitude radial velocity change will be reduced to much 
lower than 300 km/s. We suggest that the He II line comes from the X-ray photoionized 
nebula/secondary instead of the accretion disk. The line shift may be caused by noise 
given the low spectral resolution and low signal-to-noise ratios of the observations. Further 
observations that sample different orbital phases, with higher spectral resolution and higher 
signal-to-noise ratio, will help to interpret the He II line and place tight constraints on the 
black hole mass. 
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Table 1. The HST ACS observations for NGC1313 X-2 



ID 


Filter 


ExpT 


DATE 


ACor 


ZvEGA 


Zst 


VEGAmag 


j8ola2010 


HRC/F330W 


2760 


2003-11-22 


0.420 


22.904 


23.026 


22.037±0.021 


j8ol02040 


WFC/F435W 


2520 


2003-11-22 


0.277 


25.779 


25.157 


23.470±0.017 


j8ol02030 


WFC/F555W 


1160 


2003-11-22 


0.249 


25.724 


25.672 


23.625±0.026 


j8ol02010 


WFC/F814W 


1160 


2003-11-22 


0.292 


25.501 


26.776 


23.640±0.043 


j8o!06010 


WFC/F555W 


2240 


2004-02-22 


0.249 


25.724 


25.672 


23.472±0.021 



Note. - - The columns are (1) exposure ID, (2) filter, (3) total exposure in seconds, (4) 
observation date, (5) aperture correction in magnitude, (6) zeropoint for VEGAmag, (7) 
zeropoint for STmag, and (8) VEGAmag for the counterpart. 



Table 2. The optical counterpart for NGC 1313 X-2 



Object 


R.A. 


ACIS-I 
Decl. 


error 


ACS/WFC 
R.A. Decl. 


X4 

ULX a 


03:18:25.250 
03:18:22.238 


-66:36:36.04 
-66:36:03.49 


0'!17 
O'.'Ol 


03:18:24.911 -66:36:37.09 
03:18:21.900 -66:36:04.54 


X3 
ULX 6 


03:18:29.287 


-66:35:36.53 


0'.'20 


03:18:29.015 -66:35:37.57 
03:18:21.967 -66:36:04.54 



Note. — (a) This ULX position on ACS /WFC images is derived relative to the X4 position 
on ACS/WFC images, (b) This ULX position on ACS/WFC images is derived relative to 
the X3 position on ACS/WFC images. 
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Fig. I. — The identification of the ULX and X4 on the WFC F555W image, (a) X4 was 
identified with the nucleus of a background galaxy as enclosed by a C/5 circle. The smaller 
0'.'17 circle denotes the Chandra position of X4. The line segment indicates the relative offset 
of 2'.'3 between Chandra and HST positions, (b) The ULX was identified with a unique object 
labeled as 'CI' through the identification of X4. Around 'CI", the small error circle had a 
radius of Of! 2, and the large error circle Of! 5. The line segment indicating the offset between 
Chandra and HST positions was based on the identification of X4. The cross on the inner 
edge of the C/5 circle indicates the ULX position obtained through the identification of X3 
as a saturated foreground star. The object to its north was a cosmic ray only appeared in 
the 2003 F555W observation. Both panels show 7" x 7" regions. 
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Fig. 2. — Variability between two WFC F555W observations for objects within 10" of 
NGC1313 X-2. The error bar size is equal to = \Jbm\ + bm\. The error bars are 
generally smaller for brighter objects; occasionally bright objects have large error bars be- 
cause they are in unresolved stellar fields. The counterpart for NGC1313 X-2 is highlighted 
by a large circle. The histogram is for the number of objects per 0.25 in (mi — m 2 )/crAm- 
The dotted lines denote mi — m 2 = ±4o"A m - 
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Fig. 3. — The spectral energy distributions for the ULX counterpart in comparison with stan- 
dard stars of different MK spectral types. The absolute magnitudes of these spectral types are 
taken from Schmidt-Kaler (1982). All magnitudes are converted to fluxes to enable compar- 
ison between HST and UBVRI systems. The magnitudes/fluxes in filters F330W, F435W, 
F555W, and F814W were corrected for E(B-V)=0.11 mag (A = 0.561,0.455,0.335,0.196, 
respectively), E(B-V)=0.33 mag (A = 1.684,1.366,1.005,0.587, respectively), and E(B- 
V)=0.44 (A = 2.245,1.821,1.340,0.783, respectively) assuming a Galactic extinction law. 
The overall spectral shape of the counterpart dereddened with E(B-V)=0.33 mag best 
matches the shape of an 07V star. 
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F435W-F555W F555W-F814W 

Fig. 4. — The color-magnitude diagrams for the ULX counterpart and nearby bright stars. 
Both the data and isochrones are in the HST ACS/WFC VEGAmag system. The magnitudes 
are corrected for the Galactic extinction of E(B-V)=0.11 mag. (a) M F555W vs. F435W- 
F555W. The ages for the isochrones are 10 7 , 5 x 10 7 , 2 x 10 s , and 5 x 10 s years, respectively, 
(b) M F555W vs. F555W-F814W. The ages for the isochrones are 10 7 , 5 x 10 7 , 3 x 10 8 , 10 9 , 
and 3 x 10 9 years, respectively. 
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Fig. 5. — The predicted periods for different primary masses for a Z = 0.2Z & secondary 
(8.5M , 7R & ) and for an 07V secondary (3OM , 9R & ) with solar metallicity. The periods 
are calculated using the Roche lobe sizes f(q) as tabulated by Kopal, and the analytic 
approximations by Paczynski and Eggleton. 
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Fig. 6. — The orbital velocities for the primary and secondary around the mass center in 
NGC1313 X-2. A Z = 0.2Z Q secondary (8.5 M , 7 R Q ) and an 07V secondary (30M o , 9R & ) 
are considered for the system. 



Fig. 7. — True color image for a 50" x 50" region centered at the counterpart for NGC1313 X- 
2. The image was composed with F435W,F555W,F814W images from the 2003 observation. 
The stellar field is dominated by red stars with scattered blue stars. A few tens of bright 
blue stars are clustered within an irregular shape of 10" x 5" (i.e., 180pc x 90pc) to the west 
of the counterpart. North is 7.°4 clockwise from upward vertical. 



